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Abstract 

The reaction pp — > pj(7r°7r 7r°)p s has been studied at 450 GeV/c. The 7r 7r 7r° effective mass 
spectrum shows clear 77(547) and 7r 2 (1670) signals. Branching ratios for the r?(547) and 7^(1670) 
are given as well as upper limits for the decays of the u>(782), ai(1260) and a2(1320) into 37r°. 
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In a previous analysis of the centrally produced n + ix 7r° final state clear signals of the 
77(547), w(782), oi(1260), a 2 (1320) and tt 2 (1670) were observed 0. In particular, the ai(1260) 
and a 2 (1320) were observed to decay dominantly to pre. The 7T 2 (1670) was observed to decay 
to p% and / 2 (1270)7r. In order to gain more information on the decays of these states it would 
be interesting to study the 7r 7r 7r° final state since only isospin zero isobars can contribute to 
this final state. 

This paper presents new results from the WA102 experiment on the centrally produced 
7r 7r°7r° final state in the reaction 

w^P/(AV)p s (l) 

at 450 GeV/ c. The subscripts / and s indicate the fastest and slowest particles in the laboratory, 
respectively. 

The data come from experiment WA102 which has been performed using the CERN Omega 
Spectrometer, the layout of which is described in ref. 0. 

Reaction ([]]) has been isolated from the sample of events having two outgoing charged tracks 
plus six 7s reconstructed in the GAMS-4000 calorimeter, by first imposing the following cuts 
on the components of the missing momentum: |missing P x \ < 14.0 GeV/c, |missing P y \ < 
0.20 GeV/c and |missing P z \ < 0.16 GeV/c, where the x axis is along the beam direction. 
A correlation between pulse-height and momentum obtained from a system of scintillation 
counters was used to ensure that the slow particle was a proton. 

Fig. [1| shows the two photon mass spectrum (5606 events) for 67-events when the mass of 
the other two 27-pairs lies within a band around the 7r° mass (100-170 MeV). A clear 7T° signal 
is observed with a small background. Events belonging to reaction (JI|) have been selected using 
a kinematical fit (7C fit, four-momentum conservation being used and the masses of three 7r°s 
being fixed). The major background to the 7r 7T 7r° final state comes from the decay of the 
rj' and /i(1285) to i]7t 7i°. A kinematical fit has been used to remove these events. The final 
sample consists of 3590 events and has less than 2 % contamination from the 777r°7r° final state. 
The combinatorial background is also reduced by the fact that only the combination with the 
lowest x 2 is retained. 

Fig. shows the acceptance corrected 7r°7r 7r° effective mass spectrum renormalised to the 
total number of observed events. In addition to a clear 77 (547) signal there is a broad enhance- 
ment which is probably due to the 7i 2 (1670). 

The 7T°7T 7r° mass spectrum shown in fig. ^| has been fitted with a Gaussian (a = 16MeV) 
to describe the 7/(547), a Breit-Wigner convoluted with a Gaussian (er = 32 MeV) to describe 
the 7i 2 (1670) and a background of the form a(m — m t h) b exp(— cm — dm 2 — em 3 ), where m is 
the 7T 7r 7r° mass, m t h is the threshold mass and a,b,c,d,e are fit parameters. The fit is found 
to describe the data well and yields masses for the 7/(547) and 7r 2 (1670) of: 

m(77(547)) = 545 ±0.6 ±0.5 MeV, 

m(7r 2 (1670)) = 1685 ±10 ±30 MeV 

and 

r(7r 2 (1670)) = 265 ±30 ±40 MeV. 



2 



A Dalitz plot analysis of the 7r 7r°7r° final state has been performed using Zemach tensors 
and a standard isobar model ||. The analysis has assumed a, / (980), / 2 (1270) and / (1370) 
intermediate states and that only relative angular momenta up to 2 contribute. The a stands 
for the 7T7T S-wave amplitude squared, and the parameterisation of Zou and Bugg [|J has been 
used in this analysis. 

The geometrical acceptance of the apparatus has also been evaluated over the Dalitz plot 
of the 7r°7r 7r° system in 40 MeV intervals between 0.8 and 2.0 GeV. In order to perform a 
spin parity analysis the log likelihood function, Cj = log -P? (2)5 is defined by combining the 
probabilities of all events in 40 MeV 7r 7r 7r° mass bins from 0.80 to 2.0 GeV. In order to include 
more than one wave in the fit the incoherent sum of various event fractions aj is calculated: 



£ = £io g 5> p iW + (!-£<*;) (2) 




where the term (1 — J2j a j) represents the phase space background which is a free parameter 
in each bin. The negative log likelihood function (—£.) is then minimised using MINUIT ||. 
Different combinations of waves and isobars have been tried and insignificant contributions 
have been removed from the fit. The fit generates the phase space background as that part of 
the data not associated with a given wave on a bin by bin basis and one requirement is that 
this background is a smoothly varying function that shows no residual resonance structure. 

Above 0.8 GeV only the addition of the J PC = 2~ + / 2 (1270)7r° S-wave produces a sig- 
nificant change in the likelihood. We can not exclude up to 3 % contribution of the J PC = 
2~ + i / 2 (1270)7r° D-wave in the 7r 2 (1670) region. Nor can we exclude that up to 10 % of the 
7r 2 (1670) comes from a an final state. The addition of any 1 ++ wave changes the log likelihood 
by less than 1 unit and hence is classed as insignificant. However, in order to estimate an upper 
limit on the amount of di(1260) decaying to 7r 7r 7r° we include both the l ++ a7r P-wave and 
the l ++ /o(1370)7T P-wave in the final fit. The results of the final fit are shown in fig. |3[ The 
J PC = 2 _+ / 2 (1270)7r° S-wave well describes the peak in the 7r 2 (1670) region. The 1 ++ wave is 
small and flat. 

Using the acceptance corrected number of events from fits to the 7r°7r 7r° and 7r + 7r~7r° |1J] 
mass spectra the branching ratio for the 7/(547) to 7r°7r 7r° and 7T + 7r~7r° has been calculated to 
be: 



77(547) 



7r 7r 7r° 



1.35 ±0.06 ±0.09 



77(547) -> TT + TT-TT 

which is in agreement with the PDG value [H] of 1.40 ± 0.03. The branching ratio for the 
7T2(1670) to 7r 7r 7r° and 7r + 7r^7r° has been calculated in a similar manner to be: 

tt 2 (1670) -> 7r°7r 7r° 
21 ; - 0.29 ±0.03 ±0.05 



7T 2 (1670) -> TT+TT-TT 



There is no evidence for a 7r 7r 7r° decay mode of the w(782), a 1 (1260) or a 2 (1320) therefore 
an upper limit has been calculated. The masses and widths determined from the fit to the 
7r + 7r~7r° channel have been convoluted with the experimental resolution for the 7r 7r°7r° final 
state. The number of events, N, within 90 % of the predicted resonance profile has been 
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determined and the upper limit has been calculated using 2yN . For the a; (782) 

"S - M \ < 9 x 10- (90% CL) 

CJ(782) — > 71+71 7T U 

to be compared with the PDG upper limit of 3 x 10~ 4 |J which came from one experiment. 
For the ^(1260) 

Oi(1260) -> ttVtt » . ^ _ , 
i < 8 x 10~ 3 90% CL) 



ai(1260) — > 7rT7r~7T- 

if on the other hand we used the total number of events in the 1 ++ wave we would get: 

ai (1260)^W 

«Zl(l260) -> TT+TT-TT V ' 

This imposes tight constraints on the decay of the ai(1260) to isobars which have isospin 1. In 
particular, it is in disagreement with the claimed observation of an, /o(1370)7r and /2(1270)7r 
decay modes of the ai(1260) in ref. 0. This claim was based on a Dalitz plot analysis of 
the ai(1260) observed in r decays. Combining all the 1 = decays claimed in ref. [0] we 
have calculated the number of events we would expect to observe in the 7r 7r°7r° mass spectrum 
based on the number of ai(1260) observed in the 7r + 7r~7r° final state of experiment WA102 0. 
Superimposed on the 7r 7r 7r mass spectrum in fig. |4] is the ai(1260) signal we would expect to 
observe based on this number of events. As can be seen, irrespective of any results from the 
spin analysis, the number of events in the 7r 7r 7r° spectrum is much smaller than the predicted 
signal. The fact that the r\ branching ratio we have measured is in agreement with the PDG 
value indicates that our normalisation between the 7r°7r 7r° and 7r + 7r~7r° channels is correct. 
Therefore this discrepancy in the possible 1=0 decay modes of the ai(1260) could either be due 
to an overestimate of the number of ai(1260) events in the 7r + 7r~7r° final state of the WA102 
experiment or due to an error in the spin analysis performed by CLEO in ref. J7J. 

Finally, for the a2(1320) we obtain 

a 2 (1320) -> 7T 7r 7r° , . „. _ rN 

L J < 9 x 10- 3 90% CL) 

a 2 (1320) — > 7T+7T 7T U 



In summary, a study of the centrally produced 7r 7r°7r° system shows prominent signals of 
the 77(547) and 712(1670). Branching ratios for the 77(547) and 712 (1670) are given. Upper limits 
are calculated for the cj(782), ai(1260) and <22(1320) which can be used to constrain the possible 
decays of these states to isobars with isospin zero. 
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Figures 

Figure 1: The M(77) when the other two 77 pairs lie in the ir° mass region. 

Figure 2: The 7r 7r°7r° effective mass spectrum, with fit described in the text. 

Figure 3: Results of the spin parity analysis. 

Figure 4: The 7r°7r 7r° effective mass spectrum, with superimposed the number of events ex- 
pected from the a! (1260), see text. 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 4 
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